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EFFECT OF FOREST CANOPY CLOSURE ON INCOMING SOLAR RADIANCE 


C. Lisette Dottavio 
Earth Resources Branch 
Code 923 

NASA/Goddard Space Flight Center 
Green belt, Maryland 20771 

ABSTRACT 

A previous investigation into the utility of Landsat multispectral scanner (MSS) data for detecting 
gypsy moth defoliation has met with limited success. The ability to separate healthy and moderate 
defoliation is confounded by spectral similarity and topographic effects. In order to better under- 
stand the physical processes involved in defoliat.on assessment from remotely sensed data, a field 
study was designed to investigate the effect of forest canopy closure and other environmental var- 
ables on incoming solar radiation. Diffuse radiation measurements were recorded in red, infrared, 
and middle infrared wavelengths using the Mark II Three Band Field Radiometer. Results to date 
indicate that the percent canopy closure is the single most important variable affecting incoming 
solar radiation. In the visible and near infrared regions, interaction between time of day and date 
(defined later as solar zenith angle) also affect radiometric response. Aspect has only limited in- 
fluence on radiance response. These same variables do not influence middle infrared response, 
however. Unifomity of the forest canopy appears to be more important. These results are com- 
pared to Landsat MSS classification results of gypsy moth defoliation. 
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EFFECT OF FOREST CANOPY CLOSURE ON INCOMING SOLAR RADIANCE 


INTRODUCTION 

Over the last several years, personnel at NASA’s Goddard Space Flight Center (GSFC) have been 
examining the use of Landsat Multispectral Scanner (MSS) data to detect and monitor insect defoli- 
ation of hardwood forests. Initial results indicated that defoliated areas were spectrally similar to 
selected agricultural Helds in the study area. Detecti : n of defoliation by the gypsy moth cater- 
pillar was later improved using multitemporal data sets (Williams and Stauffer, 1978). Landsat data 
collected over a given forest prior to insect infestation was classified using computer-aided analysis 
techniques to identify the extent of forest cover versus non-forest cover. A second data set, collected 
over the same area after defoliation occurred was then digitally overlaid onto the forest/non-forest 
classification map. Forested areas were isolated on the image exhibiting insect defoliation and subse- 
quent analysis was limited to those areas only. 

Several processing techniques were applied to the defoliated imagery in an attempt to discriminate 
various levels of insect damage (Williams et al., 1979). In addition to maximum likelihood proce- 
dures, vegetation indices which relate the Landsat spectral measurements to various vegetation den- 
sity indicators such as green leaf biomass were applied to the MSS data. Each procedure separated 
healthy forest from heavily defoliated stands. However, the ability to separate healthy forest from 
moderate defoliation (30-60% canopy removed) proved to be more difficult because the two cover 
types have nearly identical spectral responses (Nelson, 1981 ). This problem was further compounded 
by slope and aspect. Whereas some moderate defoliation was separable from healthy and heavily 
defoliated forest on southeastern slopes, the same level of defoliation was usually not identified on 
northwestern slopes. 

In order to better understand why certain levels of defoliation were not separable using Landsat 
data, a field study was designed to investigate the radiometric properties of forests under varying 
crown densities. The objectives of this investigation were to determine the degree of change in forest 
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canopy cover detectable by radiometric measurements and to determine what other sources of vari- 
ation contributed to the radiance measured. The experiment was carried out using a Mark II 3-band 
Radiometer developed at Goddard (Tucker et al., 1981). Tlie instrument was design^ to measure 
energy in three wavelength regions that corre^XHUl to LandMt D's Thematic Mapper band TM3 
(0.63 - 0.69 /am). TM4 (0.76 - 0.90 mui) and TM5 (1.55 - 1.75#im). 

Scn-ral investigators have used similar radiometers to examine the in situ spectral properties of agri- 
aiitural crops (Kimes, et al., i980; Tucker et al., 1980). Recently, Holben et al. (1980) used the 
radiometer to evaluate a spectral method for nondestructive leaf area index determination in a trop- 
ical I'oresf by a’cording diffuse radiance at the forest floor and within the forest canopy. This 
sampling procedure was particularly valuable for studying forest radiometric properties because the 
investigators were not hindered by the difflculties associated with collecting reflectance measure- 
ments over the forest canopy. 



COLLECTION OF FIELD DATA 

Diffuse radiance was measured under the forest canopy at 48 sample plots along Doubling Gap 
Mountain near Carlisle, PA. Measurements were made throughout the gypsy moth summer feeding 
cycle in 1980 (May, June, and July), Doubling Gap is covered by mature hardwood forest consisting 
of oak (Quercus sp), hickory {Carya sp), and maple (Acer sp). The mountain is part of the ridge and 
valley physiographic region of the Appalachians, and is subject to gypsy moth attack. 

Eight 300m x 300m blocks were randomly selected along the Mountain; four on the north facing 
slope and four on the south facing slope. Within each block, six plots, 100m x 100m, were chosen 
at random. Stand age, species composition, and basal area were recorded at each plot center. Radi- 
ance measurements were taken on cloud free, low haze days, using the Mark II 3-band radiometer. 

A 30cm X 30cm barium sulfate plate was placed on the ground at the center point of the plot and 
levelled. Holding the radiometer approximately 45 cm abo/e the plate, radiance measurements were 
talcen over 5 randomly selected points on the plate (See Figure 1). By measuring light reflected from 
the barium sulphate plate, diffuse incoming radiation was recorded in each of the three wavelength 
regions (TM3, TM4, TM5). The date and time of day were noted. These measurements were used 
in subsequent statistical analyses. 

Following the radiometer readings, a 35 mm camera with a 28mm lens was placed at the plot center 
point looking directly up at the forest c nopy Photographs were taken of the forest canopy and 
were later examined to estimate percent sky versus canopy cover. In the laboratory, a 16 dot per 
square inch grid was overlain onto each photo. Dots falling within open sky were counted and 
divided by the total number of dots covering the photo to determine percent open sky. The unifor- 
mity of the canopy (i.e„ uniformly dispersed or large holes within canopy) was also interpreted from 
the nadir viewing photos. 

After all the field data were collected and photographs interpreted, the data were analyzed with 
linear regression and analysis of covariance using the 1979 Statistical Analysis System (SAS). 
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Although heavy defoliation was expected within almost aU the sample blocks, only one received 
substantial insect damage. Consequently, the complete range of canopy closures was not sampled 
throughout this exercise. Radiometric measurement were obtained for crown closures ranging be- 
tween 60 and 97 percent and at 0 percent only. Therefore, the results of this experiment can only 
be expressed as expected tendencies which need to be further substantiated by additional field 
studies. 
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GENERAL LINEAR MODEL 


Simple linear regre&iions were used to examine the effect of canopy closure on incoming sol&r radi* 
ance. Since statistical analyses require that populations be normally distributed, an arcsine trans- 
formation was applied to the pei'cent canopy closure to normalize these data. The transformation 
was performed using the following formula: 

Pt = Arcsin 

where, 

P| = transformed percent canopy closure value 
P, = actual percent canopy closure value 

The regression lines calculated for each wavelength band are plotted in Figures 2, 3, and 4. The 
general linear model for each line is defined as: 

y^ = a + ^, + € 

where, 

y = predicted radiance value for band x 
a = intercent 
0 = coefficient 

P( = transformed percent canopy closure value 
€ = error term 

The moderately high r^ values (.71, .67, .53) are indicative of the strong correlation between canopy 
closure and incoming solar radiance. However, there still remains considerable variability not 
accounted for by the model. In fact, variability tends to increase within the mid-ranges of canopy 
closure. Although a direct correlation between the field radiance values and Landsat MSS reflec- 
tance values cannot be made, the field data analyses give several indicators which arc useful to 
Landsat based studies of insect defoliation. For example, moderate defoliation (30-60% canopy 
removed) may have such high spectral variability that the class cannot be identified, even with higher 
spectral resolution, unless other environmental variables are considered. In view of this indicator, 
analyses of covariance were used to examine the effects of other variables on incoming solar radiancx*. 
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ANALYSIS OF COVARIANCE 


Envtronmenttil variables, such as aspect, canopy uniformity, and time and date of radiance measure* 
ment, were selected as covariates in the general linear model to examine the effect of these variables 
on ncoming solar radiance. Analyses of covariance for each wavelength band was performed using 
the following model; 

yx = a + U, P, + jSjA + + ^40 + ^jC + /Jj P,*A + fiy P,*A + 0 g Pt*T + 

0g P,*D + 09 P,*C +0,0 A*T + 0„ A*D + 0,2 A*C + 0,3 T*D + 

0,4 T*C + 0,5 D«C + e 

where, 

y = predicted radiance value for band x 
a = intercept 
0 j = coefficient 

P, = transformed percent canopy value 
A = aspect 

T - time of measurement 
D = date of measurement 

C = canopy distribution (i.e., uniformly dispersed or having large holes within canopy) 
e = error term 

The results of these analyses are given in Tables 1, 2, and 3. The sums of squares listed in the tables 
are those which are calculated for the particular variable being added to the model last. The values 
do not reflect incremental sums of squares for the model. The subsequent *’R>F value is a more 
appropriate measure of the variable significance. 

In each wavelength region, no singular variable appears to influence radiance respoiise. However, the 
interactions of several variables do affect response. For example, the interactions of aspect with date 
an^ with time are particularly significant for TM3 and TM4 (red and reflective infrared wavelengths, 
respectively). Further observation shows the date * time interaction to be significant for TM3 and 
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of lesser importance for TM4. These interactive variables are indicative of the effect of incident 
solar angle on radiance - a resptmse which is also evident from Landsat MSS investigations (Justice 
ct al., 1980; Huffer et al., 1975; Holben ami Justice. 1979). Another series of statistical analyses 
using solar zenith angle in l«u of date and time will be discussed in a following section. 

Two additional observations from the analyses of covariance are worthy of comment. The percent 
of canopy closure showed no influence on radiants when used in conjunction with other environ- 
mental variables. This illustrates the ability of other variables to mask out the effect of what would 
inittaily appear to be the more dominant variable. In subsequent analyses the importance of canopy 
clii^ure will jgain become evident. 

l inallv. the response of TM5 to the general linear model is contrary to both TM3 and TM4. This 
middle infrared wavelength band does not appear to be influencec by any of the interactions noted 
im viously. nor by any singular variable. Only the date * distribution interaction is significant, indi- 
V ating that tins vvavelength region may be unique in its radiometric characteristics. 
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ANALYSIS OF COVARIANCE USING SOLAR ZENITH ANGLE 


Following the initial analyses of covariance, additional analyses were run on the radiance mearatv- 
ments replacing the date and time variables with solar zenith angle. Tables 4, 5, and 6 list the analy- 
sis results for eacii wavelength band, using the general model; 

y, - a + ^,P, + ^jA + ^jZ + ^4C t t ^^P.^Z + ^7P,*C ♦ C + 0 ^\*Z + 

0^A*C ♦ d,oZ*C t e 

where. 

y^ * predicted radiance value for band 
0i = coefficient 

P| » transformed percent canopy closure value 
A * aspect 
Z - solar zenith angle 
C = canopy distribution 

The importance of percent canopy closure, solar zenith angle, and the distribution of the forest 
canopy to radiance response in TM3 and TM4 is readily apparent from the covariance tables (Table 
4 and 5). These variables and their interactions are highly significant. The results indicate that when 
using remotely sensed data to monitor forest canopy condition, solar zenith angle should be ac- 
counted for to make optimum use of radiance (or reflectance) measua'ments. Interestingly, the 
effect of aspect, seen previously, is not evident when zenith angle is accounted for. 

In the middle infrared wavelength region (TM5), solar zenith angle does not anpear to influence in- 
coming solar radiance (See Table 6). in contrast, the distribution of the forest canopy is the most 
significant variable followed by the Arcsin * Aspect and Arcsin * Distribution interactions. The rel- 
atively low r' value of 0.59 is indicative of the absence of relevant variables which account for much 
of the variability in radiance response for middle infrared wavelengths. This again points to the 
uniqueness of this spectral region. 
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CONCLUSIONS 


The technique employed in this study enabled the investigator to examine the radiometric properties 
of several different fo'-wSt conditions. Diffuse solar reliance mea»ired at the forest floor revealed 
responses to forest cover that were comparable to results obtained from previous Landsat studies. 
Several conclusions were reached, based on this study: 

1 . Percent canopy closure is the single most important variable affecting radiance within a 
forest stand. 

2. Solar zenith angle, expressed as a date * time interaction or as incident angle, influences 
radiance response in visible (red) and near infrared wavelength regions. 

3. A topographic effect, attributed to aspect, alone, in previous gypsy moth studies, was not 
evident in this study. 

4. Percent canopy closure and solar zenith angle did not appear to have significant influence 
on the middle infrared wavelength region examined. The uniformity of the forest canopy 
appeared most important. 

The conclusions reached through the results of this project are in need of further research and evalu- 
ation. Anticipated research efforts will involve a summer 1981 Held radiometer study as well as a 
spectral reflectance investigation to begin in the spring of 1981. 
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Table 1 

Analysis of Covariance for TM3 (r^ » 0.92) 


S0UT(» 

Degrees of 
Freedom 

Type IV 
Sum of Squares 

F Value 

PR>F 

Arcsin of Percent Canopy 

1 

0.63 

1.19 

0.2798 

Aspect 

1 

0.23 

0.44 

0.5091 

Date 

2 

0.28 

0.26 

0.7685 

Time 

15 

9.13 

1.15 

0.3305 

Canopy Distribution 

1 

1.06 

2.02 

0.1603 

Arcsin • Aspect 

1 

0.49 

0.92 

0.3397 

Arcsin • Date 

2 

0.16 

0.15 

0.8579 

Arcsin * Time 

15 

7.22 

1.05 

0.4153 

Arcsin • Distribution 

1 

1.40 

2.65 

0.1085 

Aspect • Date 

2 

7.51 

7.11 

0.00 16-'- 

Aspect * Time 

14 

10.61 

1.83 

0.0666+ 

Aspect * Distribution 

1 

1.20 

2.26 

0.1371 

Date * Time 

17 

15.35 

1.82 

0.0471+ 

Date * Distribution 


1.36 

1.28 

0.2840 

Time * Distribution 

13 

6.03 

0.88 

0.5799 


♦SignificMit Effect 
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Table 2 

Analysis of Covarian<» for TM4 (r^ ■ 0.9 1 ) 


Source 

Degrees of 
Freedom 

Type IV 
Sum of Squares 

F Value 

PR>F 

Arcsin of Percent Canopy 

1 

1.94 

0.82 

0.3675 

Aspect 

1 

1.54 

0.66 

0.4209 

Date 

1 

'y 

dm 

1.94 

0.41 

0.6637 

i 

1 Time 

15 

37.48 

1.06 

0.4069 

; Canopy Distribution 

1 

3.34 

1.42 

0.2379 

Arcsin • Aspect 

1 

2.95 

1.25 

0.2669 

1 Arcsin * Date 

! 

'y 

dm 

1.27 

0,27 

0.7644 

1 

Arcsin * Time 

15 

31.07 

1.02 

0.4466 

! Arcsin * Distribution 

1 

4.43 

1.88 

0.1747 

Aspect • Date 


29.16 

6.20 

0.0034+ 

! Aspect * Time 

14 

43.61 

1.69 

0.0954+ 

Aspect * Distribution 

1 

8.78 

3.73 

0.0575+ 

j Date • Time 

17 

58.93 

1.57 

0.1032 

' Date Distribution 

! 

T 

5.09 

1.08 

0.3443 

j rime • Distribution 

>3 

31.83 

1.04 

0.4246 


♦SiniitU'.inr 1 tTivt 
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Table 3 

Analysis of Covariance for TM5 (r^ » 0.82) 


Source 

Degrees of 
Freedom 

Type IV 
Sum of Squares 

F Value 

PR.'>F 

Arcsin of Percent Canopy 

1 

0.21 

0.85 

0.3609 

Aspect 

1 

0.18 

0.72 

0.3989 

Date 

2 

0.17 

0.35 

0.7072 

Time 

15 

2.79 

0.76 

0.7129 

Canopy Distribution 

1 

0.20 

0.83 

0.2785 

Arcsin • Aspect 

1 

0.29 

1.19 

0.3642 

Arcsin * Date 

2 

0.19 

0.40 

0.6743 

Arcsin • Time 

15 

2.12 

0.67 

0.7856 

Arcsin • Distribution 

1 

0.30 

1.24 

0.2690 

Aspect • Date 

2 

0.45 

0.93 

0.4008 

Aspect * Time 

.14 

2.28 

0.85 

0.5915 

Aspect • Distribution 

1 

0.02 

0.10 

0.7562 

Date * Time 

17 

4.54 

1.16 

0.3188 

Date * Distribution 

2 

1.96 

4.04 

0.0221 + 

Time * Distribution 

13 

3.98 

1.25 

0.2623 


^Significant Effect 
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Table 4 

Analysis of Covariance for TM3 uaing aolar zenith an^e in place of date and time 

(r3* 0.80) 


Source 

Degrees of 
Freedom 

Type IV 
Sum of Squares 

F Value 

PR>F 

Arcsin of Percent Canopy 

1 

8.58 

14.54 

0.0002+ 

Zenith Angle 

1 

8.30 

14.07 

0.0003+ 

Aspect 

1 

0.17 

0.28 

0.5972 

Canopy Distribution 

1 

6.36 

10.77 

0.0013+ 

Arcsin • Angle 

I 

2.99 

5.07 

0.0258+ 

Arcsin * Aspect 

1 

0.58 

0.99 

0.3216 

Arcsin "■ Distribution 

1 

4.20 

7.11 

0.0085+ 

Angie * Aspect 

1 

1.50 

2.54 

0.1132 

Angle * Distribution 

1 

3.42 

5.79 

0.0174+ 

Aspect * Distribution 

1 

0.49 

0.83 

0.3643 


♦SiKnifKant I ffcct 
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Table S 

AnalyiUs of Covariance for TM4 imng solar zenith angle in place of date and time 

(r* * 0.77) 


Source 

Degrees of 
Freedom 

Type IV 
Sum of Squares 

F Value 

PR>F 

Arcsin of Percent Canopy 

1 

20.11 

7.45 

0.0071+ 

Zenith Angle 

1 

19.33 

7.16 

0.0083+ 

Aspect 

1 

0.23 

0.09 

0.7710 

Canopy Distribution 

1 

26.35 

9.76 

0.0022+ 

Arcsin • Angle 

1 

2.08 

0.77 

0.3814 

Arcsin • Aspect 

1 

5.07 

1.88 

0.1726 

Arcsin * Distribution 

1 

16.19 

6.00 

0.0155+ 

Angle * Aspect 


4.72 

1.75 

0.1880 

Angle * Distribution 


18.85 

6.98 

0.0091 + 

Aspect • Distribution 

- - - ■ -1 

1 



4.03 

1.49 

0.2239 


^Siftnificant KfTcvt 











Table 6 

Analysis of Covariance for TMS using solar zenith angle in place of date and time 

(r2 - 0.59) 


Source 

Dt^reesof 

Freedom 

Type IV 
Sum of Squares 

F Value 

PR>F 

Arcsin of Percent Canopy 

1 

0.52 

1.98 

0.1617 

Zenith Angle 

1 

0.13 

0.50 

0.4824 

Aspect 


0.18 

0.71 

0.4014 

Canopy Distribution 

1 

1.21 

4.63 

0.0330*^ 

Arcsin * Angle 


0.001 

0.00 

0.9465 

Arcsin * Aspect 

1 

1.18 

4.52 

0.0352+ 

Arcsin * Distribution 

1 

0.93 

3.56 

0.0613+ 

Angle * Aspect 

1 

0.07 

0.26 

0.6076 

Angle * Distribution 

1 

0.44 

1.68 

0.1967 

Aspect * Distribution 

1 

0.01 

0.05 

0.8276 


+Significant Mlctt 
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1 ij-Miri.' 


I. DilTuso radiance was measured at the forest lloor usinj! ilie Mark II 3-Hand Kadiometer 

lield o»er a hanum >ulfate plate. 
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ARCSIN . PtRCtNT CANOPV COVIR 
'Rcdxtnt' 

I’ipure 2 i’ridicted radiance cufM’ lor canopy closures jirealer 
than 60 percent (0.63 - 0.6‘>/jnn. 
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Figure 3. Predic ( id radiance curve for canopy closures greater 
than 60 percent (0.76 - 0.90 Mm). 
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Figure 4. Predicted radiance curve for canopy closures greater 
than 60 percent (1.55 - 1.75 Mm). 
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